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The  local  structure  of  Ko.sNao.sNbOs  is  investigated  using  first-principles  methods  with  an  optimized  special 
quasirandom  structure  (SQS).  Through  a  comparison  of  the  computed  pair  distribution  functions  with  those 
from  neutron  powder  diffraction  data,  the  SQS  approach  demonstrates  its  ability  to  accurately  capture  the  local 
structure  patterns  derived  from  the  random  distribution  of  K  and  Na  on  the  perovskite  A -site.  Using  these 
structures,  local  variations  in  Na-0  interactions  are  suggested  to  be  the  driving  force  behind  the  R3c  to  Pm  phase 
transition.  A  comparison  between  the  SQS  and  a  rocksalt  structure  shows  the  inability  of  the  latter  to  account  for 
the  local  variability  present  in  a  random  solid  solution.  As  such,  the  predictive  nature  of  the  SQS  demonstrated 


here  suggests  that  this  approach  may  provide  insight  i 
oxide  alloys  or  solid  solutions. 
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I.  INTRODUCTION 

A  solid  solution  of  KNbOs  (KNO)  with  NaNbOs  (NNO) 
was  recently  investigated  as  an  environmentally  friendly 
alternative  to  Pb  containing  Pb(Zr,Ti)03  (PZT)  [1-3].  Similar 
to  PZT,  KNO-NNO  shows  a  significant  phase  transition  at 
a  composition  near  (Ko.5Nao.5)Nb03  (KNN)  [4].  While  the 
change  in  phase  is  not  mediated  by  a  monoclinic  phase  and 
the  piezoelectric  response  is  still  relatively  low,  understanding 
the  nature  of  the  electromechanical  properties  of  this  material 
near  this  phase  boundary  remains  important  for  understanding 
this  crystal  family. 

Given  the  relative  chemical  similarities  between  K  and  Na, 
it  is  anticipated  that  these  cations  will  be  distributed  randomly 
on  the  perovskite  A-site.  Such  random  alloys  represent  a  major 
challenge  for  first-principles  calculations  due  to  the  limited 
size  of  the  computational  supercells  that  can  be  simulated. 
Typically,  the  structure  and  properties  of  an  alloy  are  studied 
using  small  supercells,  approximating  the  alloy  compositions 
with  highly  ordered  structures  [5-7].  For  compounds  with 
simple  ratios,  like  KNN,  a  rocksalt  structure  (RSS)  is  often 
used  [8].  The  RSS  necessarily  has  short-ranged  structural 
correlations  that  deviate  significantly  from  those  anticipated 
for  a  random  structure,  and  these  may  have  consequences  for 
predicted  macroscopic  properties. 

Although  it  is  impossible  to  represent  a  truly  random 
structure  using  a  system  with  periodic  boundary  conditions, 
there  exist  several  approaches  to  approximate  the  structure 
and  behavior  of  random  alloys.  Mean-field  methods,  such  as 
the  virtual  crystal  approximation,  do  reproduce  the  properties 
in  systems  in  which  an  averaging  approach  can  be  employed 
[9,10].  However,  it  is  well  known  that  variations  in  short- 
ranged  chemical  ordering,  which  are  not  expressed  by  these 
chemical  averaging  methods,  play  an  important  role  in  defining 
the  behavior  and  properties  of  complex  oxide  ferroelectric s 
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understanding  the  properties  of  a  wide  range  of  bulk 

number(s):  77.84.-s,  61.05.fm,  61.66.Dk,  71.15.Mb 

[11-13].  Alternatively,  interatomic  potentials  that  can  model 
large  unit  cells  can  account  for  the  effects  of  local  variations 
in  bonding.  These  methods  have  been  demonstrated  to  be 
successful  at  representing  the  physical  properties  of  the  ABO^ 
perovskite  compounds  PZT  [14]  and  BaTi03-SrTi03  [15,16]. 
Unfortunately,  these  approaches  require  extensive  effort  to 
parametrize  the  potentials  and  necessitate  large  experimental 
or  computational  datasets. 

In  this  paper,  first-principles  density  functional  theory 
is  used  with  the  special  quasirandom  structure  (SQS)  [17] 
approach  to  explore  the  properties  of  the  complex  random 
alloy  KNN.  A  SQS  is  designed  by  selectively  populating  a 
supercell  such  that  the  short-ranged,  geometric  correlations 
approximate  that  of  the  random  alloy.  Although  this  approach 
has  been  demonstrated  to  accurately  reproduce  the  properties 
of  metals  [18-20]  and  semiconductors  [21,22],  it  has  only 
received  minimal  attention  by  the  oxide  community  [23]. 
The  predicted  results  from  these  computations  are  compared 
to  structures  measured  from  powered  neutron  diffraction 
experiments.  This  method  gives  pair  distribution  functions 
(PDFs)  that  are  in  excellent  agreement  with  experimental  low- 
temperature  R3c  and  room- temperature  Pm  phases.  This  is  in 
stark  contrast  with  the  ground- state  structures  obtained  using 
a  typical  RSS.  Interestingly,  the  RSS  produces  a  much  broader 
distribution  of  cation-0  and  0-0  displacement  patterns.  The 
SQS,  however,  shows  that  the  polarization  direction  and  hence 
the  macroscopic  phase/symmetry  are  correlated  most  strongly 
to  local,  nearest-neighbor,  Na-0  interactions  that  change 
most  significantly  between  the  ground- state  R3c  and  room- 
temperature  Pm  phases.  This  indicates  that  Na  interactions  are 
the  driving  force  for  the  observed  electromechanical  responses 
in  these  materials. 

II.  METHODS 
A.  Theoretical  methods 

The  structure  and  properties  of  the  SQS  and  RSS  su¬ 
percells  were  calculated  using  density  functional  theory  as 
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TABLE  L  The  lattice  vectors  and  atomic  positions  of  the  SQS  in 
both  absolute  coordinates  (scaled  by  a)  and  relative  coordinates. 


Lattice  vectors 

a 

b 

c 

(-1,1,1) 

(1-1,1) 

(2,2,0) 

Atomic  positions 

Absolute 

Relative 

K  atoms 

(2,1,1) 

(0.25,0.75,0.75) 

(1,2,1) 

(0.75,0.25,0.75) 

(0,0,1) 

(0.5, 0.5,0) 

(1,1,1) 

(0.5, 0.5, 0.5) 

Na  atoms 

(1,0,1) 

(0.25,0.75,0.25) 

(0,1,1) 

(0.75,0.25,0.25) 

(0,0,0) 

(0,0,0) 

(1,1,0) 

(0,0, 0.5) 

implemented  in  the  QUANTUM  ESPRESSO  software  package 
[24].  Ultrasoft  pseudopotentials  were  used  in  place  of  the 
all-electron  ion  potentials  [25],  and  the  exchange  correlation 
energy  was  approximated  using  the  local-density  approxi¬ 
mation  (LDA)  [26].  The  Brillouin  zone  was  sampled  using 
a  Monkhorst-Pack  [27]  6x6x6  k-point  mesh,  and  the 
wave  function  was  represented  as  a  plane-wave  expansion 
that  was  truncated  at  an  energy  of  100  Ry.  Using  these 
approximations,  the  atomic  structures  were  optimized  until 
the  residual  forces  on  the  atoms  were  less  than  5  meV/A. 
The  cubic  lattice  parameters  for  KNbOs  and  NaNbOs  were 
calculated  to  be  3.95  and  3.91  A,  which  agrees  well  with 
experiment  [28,29].  Polarization  and  Born  effective  charges 
were  calculated  from  density  functional  perturbation  theory 
[30]  and  the  Berry-phase  technique  [31]. 

The  Alloy  Theoretic  Automated  Toolkit  (atat)  [32]  was 
used  to  produce  SQSs  with  a  1  :  1  K:Na  ratio.  The  SQSs 
generated  had  eight  cubic  unit  cells  with  lattice  parameter,  a. 
The  optimal  SQS  was  selected  as  the  one  with  the  minimum 
root-mean-square  deviation  from  the  random  pairwise  corre¬ 
lations  out  to  a  distance  of  6a.  It  should  be  noted  that  although 
only  one  SQS  was  tested  here,  a  total  of  five  SQSs  were  found 
to  have  correlations  that  identically  match  the  correlations 
for  a  random  alloy  out  to  third  nearest  neighbors,  a^/3.  The 
selected  SQS  lattice  vectors  and  unit  cell  occupations  are  listed 
in  Table  I  and  depicted  in  Fig.  1 .  The  SQS  is  a  superlattice  in 
the  [110]  direction  with  a  periodicity  of  2\f2a.  The  stacking 


PIG.  1.  (Color  online)  The  atomic  supercells  used  in  these  cal¬ 
culations  to  represent  (Aq^Aq^)  5O3  perovskite  random  alloys.  The 
SQS  is  shown  on  the  left  and  the  RSS  on  the  right. 


sequence  is  . . .  ,NaNa,NaK,  KK,  KNa,. . . ,  where  NaNa  refers 
to  planes  of  all  Na  and  NaK  refers  to  planes  containing 
alternating  columns  of  Na  and  K  sites  aligned  in  the  [ill] 
direction. 

Both  the  R3c  and  Pm  phases  were  examined  [33].  In 
both  cases,  the  atomic  structures  were  initialized  according 
to  their  experimentally  reported  Glazer  tilt  patterns: 
for  R3c  and  for  Pm  [34].  The  R3c  atomic  structure 

was  allowed  to  fully  relax.  The  Pm  phase  was  constrained 
to  maintain  its  symmetry,  which  limited  the  atoms  to  moving 
only  in  the  ac  plane. 

B.  Experimental  methods 

For  comparison,  a  set  of  KNN  specimens  was  prepared  by 
solid  solution  reaction  methods  and  examined  using  neutron 
diffraction.  The  starting  materials,  Na2C03  (^99.9%,  Fisher), 
K2CO3  (^99.0%,  Alfa  Aesar),  andNb205  (^99.99%,  Sigma- 
Aldrich),  were  stoichiometrically  mixed  in  ethanol,  and  the 
resulting  composition  was  then  vibratory  milled  with  zirconia 
milling  media  for  12  h.  After  drying  for  20  h,  the  powders  were 
calcined  at  950  °C  for  5  h.  The  composition  and  phase  purity 
of  the  specimen  were  verified  using  a  Siemens  Model-D500 
x-ray  diffractometer.  The  resulting  powder  x-ray  diffraction 
pattern  for  the  KNN  is  shown  in  Fig.  2. 

High-  Q  neutron  powder  diffraction  data  were  collected  on 
the  NPDF  instrument  at  the  Los  Alamos  National  Laboratory 
Lujan  Neutron  Scattering  Center  [35].  Approximately  10  g 
of  KNN  powder  was  loaded  into  a  vanadium  can,  and 
measurements  were  performed  at  15  and  300  K  to  allow  for 
examination  of  both  the  rhombohedral  and  monoclinic  phases. 
Using  the  program  PDEGETN  [36],  standard  data  correction  and 
normalization  steps  were  performed,  and  the  experimental  pair 
distribution  functions  were  extracted.  Based  on  inspection  of 

noise  levels  in  the  data,  a  wave- vector  cutoff  gmax  of  34  A 
was  used. 

III.  RESULTS  AND  DISCUSSION 

Reduced  pair  distribution  functions  (PDFs)  from  the  ex¬ 
perimental  neutron  powder  diffraction  data  and  from  the 


(deg.) 


FIG.  2.  The  6-26  powder  x-ray  diffraction  intensities  taken  for 
the  calcined  (Ko.5Nao.5)Nb03  powder.  The  pattern  indicates  a  single 
perovskite  phase  and  is  indexed  based  on  the  primitive  cubic  cell. 
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EIG.  3.  (Color  online)  A  comparison  of  the  experimentally  deter¬ 
mined  reduced  PDE  and  the  theoretical  PDEs.  The  top  frame  shows 
the  results  for  the  rhombohedral  phase  and  the  bottom  frame  shows 
results  from  the  monoclinic  phase. 


theoretical  structures  are  compared  in  Fig.  3.  The  theoretical 
structures  are  linearly  scaled  by  1%  to  correct  for  the  well- 
known  LDA  errors  in  lattice  constants.  The  top  frame  of 
Fig.  3  shows  results  for  the  R?>c  phase.  The  peak  positions 
and  relative  intensities  of  the  SQS  are  in  excellent  agreement 
with  the  experimental  PDF  obtained  at  15  K.  In  contrast,  the 
RSS  PDF  has  spurious  peaks  and  in  several  instances  there  are 
missing  peaks.  The  results  for  the  Pm  phase,  seen  in  the  bottom 
frame  of  Fig.  3,  further  verify  the  ability  of  the  SQS  to  predict 
the  local  structure  of  the  KNN  solid  solution.  The  agreement 
between  the  SQS  and  experiment  is  again  emphasized  using 
the  refined  PDFs  as  shown  in  Fig.  4,  in  which  the  experimental 
data  are  modeled  using  the  refined  SQS  structures  and  the 
PDFGUI  software  [37].  The  refined  PDFs  clearly  demonstrate 
the  capacity  of  the  SQS  to  provide  a  structural  model  for 
the  KNN  random  alloy.  Conversely,  the  RSS  shows  obvious 


FIG.  4.  (Color  online)  Refined  PDFs  for  the  R3c  (top)  and  Pm 
(bottom)  phases. 


23456789 
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FIG.  5.  (Color  online)  Partial  pair-distribution  functions  for  all 
atomic  pairs.  The  solid  black  line  is  the  SQS  rhombohedral  phase  and 
the  dashed  red  line  is  the  RSS  rhombohedral  phase.  The  monoclinic 
phase  for  the  SQS  is  shown  as  the  dotted  green  line. 


deficiencies  that  may  be  an  indication  of  a  poor  representation 
of  the  random  alloy. 

An  analysis  of  the  SQS  and  RSS  is  performed  by  computing 
partial  PDFs,  shown  in  Fig.  5.  This  allows  for  an  understanding 
of  the  pairwise  relationships  of  these  structures,  as  well  as 
the  underlying  driving  force  behind  the  electromechanical 
responses.  Figures  5(a)-5(c)  shows  the  A-A-site  cation  in¬ 
teratomic  distances.  The  large  gaps  between  the  peaks  in  the 
RSS  data  are  due  to  the  RSS  ordering,  that  is,  there  exist 
no  nearest-neighbor  cations  of  the  same  composition,  and 
all  second-nearest  neighbors  have  the  same  composition.  The 
relative  A-A-site  cation  interatomic  distances  are  similar  in  all 
structures.  In  particular,  K-K  and  Na-Na  interactions  all  have 
narrow  distributions  indicating  strong  repulsion  between  these 
ions,  in  other  words,  collective  motion,  which  undoubtedly  sets 
the  lattice  constant  of  the  KNN  structure.  Similar  consideration 
of  the  nearest-neighbor  Nb-Nb  separations  suggests  a  lattice 
parameter  of  3.97  A  with  a  standard  deviation  of  0.015  A. 
This  is  slightly  larger  than  the  predicted  values  for  the  bulk 
cubic,  paraelectric,  NNO,  and  KNO,  and  is  a  consequence  of 
the  polarization  that  develops  in  these  unit  cells.  Nevertheless, 
in  both  the  R3c  and  the  Pm  structures  there  is  a  splitting 
of  the  nearest-neighbor  Na-K  peak.  This  is  a  result  of  larger 
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off-center  displacements  of  Na  cations  than  K.  Interestingly, 
there  is  a  shift  in  the  relative  peak  separation  between 
the  two  phases,  as  well  as  the  presence  of  a  third  peak 
in  the  Pm  phase.  This  is  likely  due  to  changes  in  the 
off-centering  of  the  two  sublattices,  which  induce  physical 
changes  in  the  direction  of  the  macroscopic  polarization  of  the 
material. 

Similarly,  sharp  regularly  spaced  peaks  are  observed  for  A- 
B  and  Nb-Nb  interactions,  as  shown  in  Figs.  5(d)-5(f),  for  both 
the  R?>c  and  Pm  phases.  Collectively,  this  is  indicative  of  a 
global  response  of  the  cations,  all  displacing  along  the  direction 
of  the  bulk  polarization.  Such  a  global  response  of  cations  to 
the  overall  polarization  has  been  previously  used  to  understand 
why  the  three  composition-dependent  phases  of  PZT  have 
very  similar  PDFs  but  very  different  macroscopic  polarizations 
[14].  In  PZT,  the  overall  macroscopic  polarization  was  found 
to  be  due  to  a  competition  between  local  distortions  due 
to  ionic  repulsion,  primarily  between  Pb  and  Zr,  and  the 
cooperative  displacement  of  the  cations  along  the  overall 
polarization  direction  [14].  Previous  work  of  Suewattana  and 
Singh  indicates  that  in  KNN,  the  5-5-site  response  is  a  result 
of  stiff  force  constants  between  the  5-site  cations  [38].  This  is 
a  direct  result  of  5-0  interactions  where  movement  of  5 -site 
cations  toward  an  O  anion  induces  the  neighboring  5 -cation 
to  displace  [39].  Unlike  PZT,  the  competition  between  cations 
is  minimized  and  only  the  collective  motion  of  the  ions  defines 
the  macroscopic  polarization.  This  could  be  a  reason  for  the 
low  piezoelectric  response  observed  in  these  materials,  since 
it  could  be  anticipated  that  the  ionic  responses  to  changes  in 
the  direction  of  the  polarization  would  be  far  less  than  those 
found  in  PZT,  where  the  changes  in  off-centering  are  much 
more  dramatic. 

Finally,  the  cation-oxygen  and  oxygen-oxygen  bonding 
interactions  are  analyzed,  as  shown  in  Figs.  5(g)  and  5(h). 
Interestingly,  partial  PDFs  are  similar  for  both  the  53c  and 
Pm  phases  for  the  interactions  between  K-0,  Nb-0,  and  0-0. 
However,  large  deviations,  especially  for  the  nearest-neighbor 
distances,  are  observed  in  the  Na-0  PDFs  between  the  53c 
and  Pm  phases.  Specifically,  the  three  nearest-neighbor  peaks 
in  the  53c  and  the  Pm  phase  are  indicative  of  different 
displacement  patterns.  In  the  53c  phase,  displacements  are  in 
the  [111]  direction,  resulting  in  three  short  bonds  (^2.65  A), 
six  medium  length  bonds  (2.65-2.9  A),  and  three  long 
bonds  (^2.9  A)  along  this  direction.  On  the  other  hand, 
in  the  Pm  structure,  the  splitting  in  Na-0  distances  is 
manifested  through  the  creation  of  one  short  bond  along 
the  [101]  direction  (^2.5  A),  four  medium  length  bonds 
(2.5-2.75  A),  and  seven  long  bonds  (^2.75  A).  Despite  these 
differences,  the  Na-0  distances  are  governed  by  the  preferred 
displacements  of  Na  toward  neighboring  cells  with  K  ions 
and  away  from  those  with  Na  cations.  In  other  words,  short 
Na-0  bonds  are  observed  in  directions  that  are  Na-O-K  and 
longer  bonds  are  in  the  direction  of  Na-O-Na  bonds.  This 
is  also  evident  in  the  splitting  of  the  Na-K  peaks.  Recent 
DFT  studies  comparing  the  cubic,  paraelectric  phases  of 
KNO  and  NNO  indicate  significant  changes  in  covalency 
and  ionic  charge  transfer  in  KNO  compared  with  NNO,  with 
K  resulting  in  more  charge  transfer  to  oxygen  ions.  This 
enhances  the  bond  hybridization,  which  fosters  off-centering 
and  induces  large  cooperative  displacements  in  nearby  cations. 


The  observed  changes  in  splitting  are  directly  linked  to  the 
desire  to  displace  in  the  direction  of  the  overall  polarization 
(or  symmetry  of  the  unit  cell).  This  cooperative-competitive 
effect  is  very  similar  to  that  observed  in  PZT,  and  tuning  this 
interaction  may  be  useful  in  enhancing  the  ferroelectric  nature 
of  KNN. 

It  is  important  to  point  out  that  the  results  for  the  53c 
obtained  using  the  SQS  method  are  in  stark  contrast  with 
those  of  the  RSS.  A  number  of  fundamental  differences  exist 
that  change  the  interpretation  of  the  origin  of  ferroelectric 
polarization  in  these  materials.  First,  for  the  Na-0  interactions 
there  are  only  two  bond  distances.  Although  the  displacements 
are  still  in  the  [111]  directions,  the  lack  of  a  third  Na-0 
bond-distance  splitting  is  indicative  of  a  smaller  magnitude  of 
off-centering  in  this  direction,  shown  in  Fig.  5(g).  Interestingly, 
there  is  now  a  more  pronounced  splitting  in  the  K-0  bond 
lengths,  as  demonstrated  in  Fig.  5(h),  which  suggests  that 
the  K  cations  now  contribute  more  to  the  ferroelectricity 
in  KNN,  that  is,  larger  off-centering,  than  indicated  by  the 
SQS.  In  addition,  although  there  is  no  change  in  the  com¬ 
puted  Nb-0  nearest-neighbor  distances,  there  is  a  significant 
broadening  of  second-nearest-neighbor  Nb-0  bonds  in  the 
RSS  versus  the  SQS,  which  is  shown  in  Fig.  5(i).  This  is 
most  likely  linked  to  the  large  variation  in  the  0-0  bond 
distances,  demonstrated  in  Fig.  5(j).  This  large  variation 
is  a  strong  indication  of  larger  octahedral  rotations.  Such 
rotations  are  usually  thought  to  counter  the  ferroelectric  modes 
in  these  materials.  Nevertheless,  these  observations  suggest 
that  by  increasing  the  chemical  complexity  of  the  supercell, 
through  a  better  approximation  for  the  randomization  of 
A -site  cations,  the  variation  in  the  bond-length  distribution 
decreases. 

Table  II  lists  the  Born  effective  charges,  Z*,  and  the 
magnitude  of  the  polarization,  5,  for  the  53c  and  Pm 
phases.  The  SQS  method  allows  for  accurate  determination 
of  the  polarization.  The  contribution  of  each  cation  to  the 
net  dipole  moment  is  characterized  by  determining  the  ions’ 
displacement  relative  to  the  center  of  the  nearest-neighbor 
oxygen  cages,  listed  in  Table  II.  On  average,  both  SQS 
and  RSS  provide  comparable  results  to  the  overall  polar¬ 
ization.  The  variability  of  ionic  displacement  due  to  the 
local  chemical  environment  is  shown  in  Fig.  6.  Cations  in 
the  SQS  undergo  larger  relative  displacements  perpendicular 
to  the  bulk  polarization  direction,  as  would  be  expected 
in  a  random  solution,  yet  these  perpendicular  components 
sufficiently  cancel  one  another  to  allow  the  greatest  contri¬ 
bution  to  be  those  components  parallel  to  the  macroscopic 
polarization.  The  RSS  shows  highly  uniform  deviation  from 


TABLE  11.  The  Born  effective  charges,  Z*,  in  units  of  electron 
charge,  angular  deviation  of  cation  displacement  away  from  the  bulk 
polarization  direction,  6,  in  degrees,  and  polarization,  5,  in  C/m^. 


Z* 

0 

P 

K 

Na 

Nb 

K 

Na 

Nb 

RSS  (53c) 

1.13 

1.14 

8.33 

7.1 

6.7 

8.4 

0.38 

SQS  (53c) 

1.15 

1.15 

8.38 

8.2 

9.3 

7.8 

0.39 

SQS  (Pm) 

1.15 

1.15 

8.33 

6.5 

9.7 

7.4 

0.40 
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EIG.  6.  (Color  online)  Distribution  of  cation  distortion  angles 
away  from  the  bulk  polarization  direction  of  supercells  for  R3c  and 
Pm  phases.  The  SQS  supercells  (broken  lines)  demonstrate  greater 
variability  in  ionic  contributions  to  overall  polarization,  while  the  RSS 
supercell  (solid  line)  is  highly  uniform  in  its  angular  distribution. 

individual  cations  due  to  the  high  symmetry  imposed  on  the 
structure.  The  SQS  is  able  to  predict  the  bulk  polarization 
while  accounting  for  local  influences  due  to  the  random 
environment. 

The  angular  momentum  and  site-projected  density  of  states 
for  the  relaxed  structures  are  shown  in  Fig.  7.  The  valence 
band  is  primarily  occupied  by  O  /7-states,  and  most  of  the  Nb 
J- states  exist  in  the  conduction  band.  Hybridization  between 
the  Nb  J-states  and  O  /7-states  is  evident  from  the  presence 
of  some  Nb  J-states  within  the  valence  band.  The  LDA  band 
gaps  for  the  ground-state  SQS  and  RSS  are  1.93  and  2.08  eV, 


Energy  (eV) 


EIG.  7.  (Color  online)  The  projected  density  of  states  for  the 
KNN  relaxed  structures.  The  origin  has  been  placed  at  the  top  of 
the  valence  band. 


respectively,  while  the  Pm  structure  has  a  band  gap  of  1.87  eV, 
in  comparison  to  an  experimental  band  gap  of  3.18  eV  [40]. 


IV.  SUMMARY  AND  CONCLUSIONS 

In  conclusion,  the  R3c  and  Pm  phases  of  the  KNN  random 
alloy  are  investigated  using  flrst-principles  calculations.  Using 
the  SQS  approach  to  construct  supercells  of  KNN,  one  obtains 
PDFs  that  are  in  excellent  agreement  with  those  obtained 
from  neutron  diffraction  experiments.  This  is  in  stark  contrast 
with  results  for  PDFs  produced  from  supercells  with  the 
RSS.  Using  the  SQS  method,  speciflc  insights  may  be  gained 
into  the  local  structure  interactions  that  control  the  physical 
properties  of  these  materials.  These  results  suggest  that  the 
ferroelectric  polarization  in  these  materials  is  governed  by  the 
cooperative-competitive  interactions  of  Na.  In  particular,  Na 
prefers  to  off-center  in  the  direction  of  K  cations  and  away 
from  Na  ions,  while  conforming  as  best  as  possible  to  the 
overall  polarization  direction.  These  interactions  seem  to  be 
fundamental  in  perovskite  oxide  alloys  and  play  an  important 
role  in  deflning  the  magnitude  and  direction  of  the  macroscopic 
polarization.  As  such,  these  results  highlight  the  effectiveness 
of  a  SQS  for  predicting  the  structures  and  properties  of  oxide 
solid  solutions  directly  from  flrst-principles  calculations  using 
relatively  small  supercells  without  resorting  to  the  use  of 
mean-fleld  methods  or  parametrized  interatomic  potentials. 
This  is  paramount  to  the  design  of  novel  materials  where 
macroscopic  properties  are  strongly  linked  to  the  local  atomic 
structure.  The  artiflcial,  high-symmetry  order  of  the  RSS  casts 
some  doubts  about  using  these  wrongly  correlated  structures  to 
predict  the  properties  of  perovskite  compounds,  although  this 
approach  is  common.  These  results  have  clear  consequences 
for  flrst-principles  discovery  of  novel  materials  with  properties 
such  as  ferroelectricity,  piezoelectricity,  ferromagnetism,  and 
thermoelectricity.  For  example,  they  suggest  that  the  response 
of  KNN  may  be  enhanced  through  modiflcations  that  increase 
the  bond  frustration  in  these  materials.  One  possibility  may 
be  through  chemical  pressure  induced  by  the  introduction  of  a 
larger  cation  such  as  Ta  [41]. 
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